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What are the constraints placed on the constitutive tensors of elastodynamics by the requirements 
that the linear elastodynamic system under consideration be both causal (effects succeed causes) and 
passive (system doesn’t produce energy)? The analogous question has been tackled in other areas 
but in the case of elastodynamics its treatment is complicated by the higher order tensorial nature 
of its constitutive relations. In this paper we clarify the effect of these constraints on highly general 
forms of the elastodynamic constitutive relations. We show that the satisfaction of passivity (and 
causality) directly requires that the hermitian parts of the transforms (Fourier and Laplace) of the 
time derivatives of the constitutive tensors be positive semi-definite. Additionally, the conditions 
require that the non-hermitian parts of the Fourier transforms of the constitutive tensors be positive 
semi-definite for positive values of frequency. When major symmetries are assumed these definiteness 
relations apply simply to the real and imaginary parts of the relevant tensors. For diagonal and 
one-dimensional problems, these positive semi-definiteness relationships reduce to simple inequality 
relations over the real and imaginary parts, as they should. Finally we extend the results to highly 
general constitutive relations which include the Willis inhomogeneous relations as a special case. 
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I. INTRODUCTION 

Various aspects of nature are modeled as cause-effect relationships between different physical processes. These 
physical processes are often functions of time and the relations between them, in some cases, can be more easily 
analyzed in the frequency or Laplace domains. A frequency dependent process is called dispersive and can be studied 
by deriving the appropriate dispersion relations of the system. Physical systems often have an inherent assumption 
of causality wherein effects are assumed not to precede causes. If the physical system is also linear and time-invariant 
then certain sum/integral rules could be derived connecting the physical quantities involved^i^. For eg. the Kramers- 
Kronig (K-K) relationships^ii are integral relationships which connect the real part of the electromagnetic index of 
refraction to its imaginary part, thus connecting dispersion and loss in the medium. Since their introduction, the K-K 
relationships have been used in the study of circuit theory^ and all forms of wave propagatioii^^— . 

The K-K relations have recently attracted interest in the area of metamaterials where the goal is to create materials 
with exotic electromagnetic, acoustic, and/or elastodynamic properties. The essential ideas emerge from early theo¬ 
retical works of Veselago^^ and more recent experimental efforts by various research groups^^^— (See^^ for a review). 
The possibility of creating materials with unprecedented material properties has led to far-reaching postulations of 
their applications, most visibly, in the area of cloakingi^r— . Since material properties can essentially be viewed as 
time-domain transfer functions which relate a cause to its effect (and, therefore, must be causal), K-K relations and 
their derivatives can be used to place some realistic constraints on the properties themselves. The K-K relations have, 
of late, been used as a tempering check on the optimism that has emerged in the area of metamaterials research. This 
causality check includes, on one end of the spectrum, placing some realistic constraints on the application potential of 
metamaterials as cloaking devices^l to, on the other end, sobering realizations that a considerable amount of meta¬ 
materials research stands on shaky foundations, often proposing materials which violate such basic ideas as causality 
and/or the second law of thermodynamics^. This has led to a number of researchers advocating a need for improved 
models for metamaterials^^—. 

Closely connected to the idea of causality is the concept of passivity which refers to the assumption that the physical 
process under consideration cannot produce energy^^. In fact, if the physical process (cause-effect relationship) can 
be expressed in a convolution form in the time domain then its satisfaction of the passivity requirement automatically 
means that it satisfies causality as well^^. A physical process can, in turn, be expressed in the convolution form if 
it satisfies certain conditions such as linearity and time-invariance. It becomes interesting, therefore, to understand 
what constraints are placed upon a linear time-invariant cause-effect relationship (constitutive relationship) in electro¬ 
magnetic, acoustic, and/or elastodynamic areas by the requirement of passivity. Such knowledge can be used to place 
constraints on and understand the limitations of various metamaterial models which are used in these areas to arrive 
at such relationships. Considerable research in this direction has already taken place in the field of electromagnetics 
where it is clear that passivity demands that the imaginary parts of the diagonal values of the Fourier transform of 
e, /X be non-negative for all positive values of frequency^^— (fields assumed to depend upon However, it is not 

clear, to the author’s knowledge, what should be the equivalent constraints in elastodynamics for the most general 
constitutive cases. The case of 1-D longitudinal or shear wave propagation in an elastic medium is equivalent, in form, 
to the electromagnetic case. As such, it immediately follows that passivity should require that the corresponding 
1-D material properties (modulus and density) should behave analogously to the e, pL. However, in 2- and 3-D, the 
elastodynamic constitutive tensor cannot, in general, be diagonalized. Moreover, recent advancements^Sf^^— suggest 
that the Willis constitutive relation^, which is a coupled form of constitutive relation, is more appropriate for the 
description of inhomogeneous elastodynamics and, therefore, of elastodynamic metamaterials. It is not clear what the 
constraints of passivity are on such highly general elastodynamic constitutive forms. 

In this paper we study the constraints which passivity places on highly general forms of the elastodynamic 
constitutive relations. We use the passivity condition which is equivalent to the statement of passivity used in 
electromagnetism^^ and circuit theory^! and which is elaborated in subsequent sections. We also present our analysis 
within the context of distributions which is the proper space within which to describe the transfer functions of passive 
systems. Furthermore, treating the constitutive tensors in the space of distributions ensures that the analysis applies 
to the metamaterial cases of most interest and also to the static case (elastic case). 


II. BACKGROUND 

Physical processes in the real world are often described as an interplay between physical variables and fields which 
are dependent upon time. The relationships between the physical variables can be modeled as input-output relations 
where a time dependent variable v{t) is produced from another time dependent variable u{t) through some rule 
IH, v{t) = iyiu{t). Although physical variables often satisfy certain continuity and differentiability conditions, it is 
desirable to consider them as generalized functions for broader applicability. We identify four spaces at this point. 
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Space X> is the space of all complex valued functions which are infinitely smooth and with compact support. It 
is a subset of space S which consists of all infinitely smooth and complex-valued functions called functions of 
rapid descent, such that they and all their derivatives decrease to zero faster than every power of l/|t| as |t| —)• oo. 
Space S of distributions of slow growth is the space of continuous linear functionals on S. Space V is the space of 
continuous linear functionals on T). It can be shown that T) C S G S GT) . Input-output relations can be completely 
arbitrary but they reduce to a particularly simple form if the properties of single-valuedness, linearity, continuity, 
and time-translational invariance are assumed to hold for the operator These conditions are generally true for 
physical processes and under these conditions the operator 91 reduces to a convolution operation: 

v{t) = R*u = / dr R{t — t)u{t) (1) 

Jr 


where the last equality only holds if R, u are locally integrable distributions whose supports satisfy certain boundedness 
properties (either R,u have bounded supports, or both R,u are either bounded on the left or on the right). The 
operator R is causal if it is not supported on t < 0. The final property of passivity can be stated by defining the 
energy of the system. If the power absorbed by the system at time s is given by Re v*{s)u{s) where * denotes the 
complex conjugate, then define the energy absorbed by the system up to time t as: 

f(t)=Re f dsv*{s)u{s) (2) 

J —OO 

The operator R is considered passive if £{t) >0 for all f S K. For operators in convolution passivity implies causality 
and it also implies that R € S . Therefore, Eq. @ is well defined for at least for all u G 5. For distributions in S 
the Fourier transform, denoted by must satisfy {J-f, </>) = (/, Tcj)) for / G 5 , (/> G 5 where (/, 0) is the value in C 
that / assigns to (p through the operation f{t)(j){t)dt. Furthermore it can be shown that the Fourier transform 

of a distribution which is in S is itself in S . For distributions in S the Fourier transform is defined in the usual way 
(w G K): 

(j) = (j)(w) = [ dt Pit) =-^ [ dw (3) 

Jr Stt Jr 

In addition to R being a distribution of slow growth, causality implies that its support is in [0,oo). For such a case 
its Laplace transform is given by i? = (i?, e~^*) and it can be shown that passivity directly implies that Re i? > 0 for 
all Re z > 0 (passivity constraint) and that R is analytic there^. It is important to point out here that if R is causal 
and is in S then its Laplace transform can be derived from its Fourier transform through its analytic extension in the 
right half of the complex plane. This is the point of departure for deriving K-K relations in electromagnetism when 
the constitutive tensors can be expressed diagonally. 

For transfer functions of higher order and more complexity which form the object of study of this paper, we 
need to dehne some additional spaces. We will use bold symbols to denote tensors whose elements are distributions. 
If f(f) is a tensor of distributions, is the matrix of complex numbers obtained by replacing each element 

of f(t) by the number that this element assigns to the testing function p{t). We will use additional subscripts with 
the spaces already defined above to denote the space in which all tensors of the relevant rank and distribution lie. 
For e.g. Rnxnxnxn is the space of all fourth order tensors whose elements are distributions in V etc. For a tensor 
of even rank f we also define the operations f^ which denotes a transpose over the major symmetry and ff which 
denotes a transpose over the major symmetry followed by conjugation. Now a single-valued, linear, time-invariant, 
and continuous input output relation can be written in the convolution form: 


v(t) = R * u = 


dr R(f — t)u(t) 


( 4 ) 


where appropriate boundedness of R, u are assumed and where v is a tensorial quantity derived from u through the 
linear operator R. Total energy absorbed up to time t is given by: 

f(t) = Re f ds v'l'(s)u(s) (5) 

J —OO 

The operator R is considered passive if £{t) >0 for alH G M and to make sure that the integral above exists, we will 
restrict the elements of the input u to be in If the operator R which is in the convolution form is passive then 
it can also be shown to be causal^^ and, furthermore, its elements are in S . Moreover, certain important properties 
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of the Laplace and Fourier transform of R follow. Its Laplace transform is given by R = (R, e and it is analytic 
(its elements are analytic) in the open right half of the complex plane. It can be derived from its Fourier transform 


through analytic continuation in its region of convergence. 


And finally, its hermitian part (R^ = .5 


R + Rt 


IS 


positive semi-definite for all Re z > 0. Note that these results follow from some fairly unrestrictive constraints on the 
input field and transfer function which are easily satisfied in elastodynamics (and electromagnetism). Our effort here 
is to apply and extend these results to the elastodynamic case. 


III. CAUSALITY AND PASSIVITY IN ELASTODYNAMICS 

We begin by considering a volume O within which the pointwise elastodynamic equation of motion and kinematic 
relations are specified: 


^ij,j 4“ fi —Pii ^ij — (6) 

where cr, e, p, u and f are the space and time dependent stress tensor, strain tensor, momentum vector, displacement 
vector, and body force vector respectively. These relations need to be supplied with appropriate constitutive relations 
which relate the various field variables to each other. For the current discussion we consider stress and velocity to 
be independent fields (input fields) which lead to the emergence of strain and momentum fields (output/dependent 
fields) respectively. The relationships are expressed in terms of general constitutive operators whose properties need 
to be determined based upon the various subsequent assumptions about the system: 

e(x,t) = 2)(cr(x,t)), p(x,t) =‘P(u(x,t)) (7) 


Now we assume that the operators satisfy the conditions of single-valuedness, linearity, time-invariance, and conti¬ 
nuity and, furthermore, that the stress and velocity fields along with the operators satisfy appropriate boundedness 
conditions referred to in the last section. We also assume that the operator is real valued, i.e. it assigns a real output 
field to a real input field (but can assign a complex output field to a complex input field). Under these assumptions 
the constitutive relations of Eq. ([7]) can be specialized to the following form: 


e(x, t)=D*cr= / ds D(x, t — s) : (t(x, s) = / ds D(x, s) : cr(x, t — s) 
J —oo j —oo 

/ GO poo 

ds p(x, t — s) • u(x, s) = / ds p(x, s) • u(x, t — s) 

-OO J —GO 


( 8 ) 


where the components of the constitutive tensors D, p are real valued distributions. D is a fourth order tensor field 
in and possesses the usual minor symmetries that are associated with stiffness and compliance tensors, p 

is a second order tensor field in To make sure that energy, as defined later on, exists, we will take cr to be in 

Rnxn and ii to be in Rnxi- No major symmetries are assumed at this point. We also note a further result which will 
be used later. For convolutions of distributions as appearing above we note that e = D * cr and p = p * ii. 

Causality: Causality refers to the requirement that an effect cannot precede its cause. With reference to the 
constitutive relations in Eq. (jS)), it implies that the value of the strain and momentum fields at time to can only 
depend upon the values, respectively, of the stress and velocity fields at times prior to and including tg- A necessary 
and sufficient condition for a system to be causal is that its unit response function (constitutive operator in the present 
case) vanishes for t < 0. Specifically, causality implies the following for the constitutive tensors: 


D(x, t) = 0, p(x, t) = 0 V< < 0 


(9) 


Passivity: Passivity refers to the requirement that the system cannot generate energy. For the elastodynamic case 
the total energy at any time t contained in U comprises of the elastic energy contribution and the kinetic energy 
contribution: 


f(t) = Re - / dx [(t(x, t) : e*(x, t)-f p(x, t) • u*(x,t)] 
2 Jn 


( 10 ) 


The total energy can, therefore, be given by the following integral: 


S{t)= J J dx ^ [(t(x,s) : e*(x, s) -bp(x, s) ■ u*(x, s)] 


( 11 ) 
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Passivity requires that the total energy absorbed by the system be non-negative at all times: 

£{t) > 0 Vt 


(12) 


To calculate the total energy absorbed by the system in Eq. m, we note that the time derivative of £{t) should 
equal the power input from the tractions, t(x, <), which are acting on d£l and the body forces, f, which are acting in 
This power input is given by: 

P(t) = Re f dxti{x,t)u*{x,t) + f dx fi{yL,t)u*{yi,t) = Re f dx [{aijU*)j + fiii*] (13) 

.J on J Q. . J Q. 

where the dependence on space and time is suppresed. The above is achieved through the application of the Gauss 
theorem and the relation ti = Uijtij. By decomposing -I- where e and w are the symmetric and 

antisymmetric parts of respectively, and noting that the inner product of a symmetric tensor and an antisymmetric 
tensor goes to zero, we have: 

P{t) = Re [ dx + f,u*] 

JQ 

Using the equation of motion (Eq. 15]) and rearranging, we have: 

P(t) = Re dx + Piii*] 

Jn 

The above calculated power should equal the rate of change of the total energy stored in Q: 

d£{t) 


(14) 


(15) 


dt 


= P(t) 


Eqs. (I11I12I16I) together give: 


£{t) = [ ds P(s) = Re [ ds [ dx [crij(x,s)e* (x,s)-|-pi(x,s)M*(x,s)l > 0 

J— oc J— oc Jq 


(16) 


(17) 


Since £l is arbitrary, the above would be satisfied only if the inequality holds at each point in space. Going forward 
we, therefore, understand the passivity statement to be the following: 


Re 


ds [cry(s)eb(s)-bp,(s)u*(s)] >0 


(18) 


where the variables are being evaluated at an arbitrary but same location x G fl. Since the constitutive tensors in 
Eq. (|8]) are real we note that e* = D * cr* and p = jj* u. Using these relations we have: 


Re 


ds 


(Tij{s) I dv bijki{v)ali{a - v)+u*{s) I dv pij{v)uj{s - v) 


> 0 


(19) 


We note that the time derivatives of the constitutive tensors should themselves be real and causal and, therefore, 
must satisfy D(x, t) = p(x, t) = 0 for all t < 0. Now we employ the distributional Laplace transform. This is done by 
first choosing <t(s) = ( T ( f >*( s ) and u(s) = U7(s) where (t,u are constant tensors and (j){s),"f{s) are in T>: 


Re 


ds 


(Tij(j)*{s) dvbijkiiv)ali(l){s-v)+u*'y*{s) dv pij{v)uj'y{s - v) 


> 0 


( 20 ) 


and then letting be equal to for —oo < s < a where t < a < oo and z G C. With the requirement of 

causality on D,p the above inequality reduces to the following: 


Re 


or 


o'ilCTfcz / dv Dijki(v)e ''''+u*Uj / dv pij{v)e 
Jo Jo 


POO POO 

(Jijali / dv bijki{v)e~^^ +u*Uj / dv pij(v)e~ 


> 0 


> 0 


( 21 ) 


Re 


( 22 ) 
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Since the constitutive tensors are distributions of slow growth, the integral quantities are identified as the distributional 
Laplace transforms of the respective constitutive tensors so that the above relations can be expressed in the following 
condensed form; 


Re 


(T 


t> : cr* +u ~p u* 


> 0 


(23) 


We now decompose the tensors D and p into their hermitian and non-hermitian parts: 

- + i . 

2 L J 2 1 


D = +D"^; 




'p=''ph + ''pnh. 


h^ = l[h + h^ 


■jUh ^ 


1 


D-Dt 
p — p^ 


2 L' 


(24) 


and note that only hermitian parts of the tensors contribute to the real part of Eq. (1^51) . Furthermore, since the 
tensors cr, li in Eq. (1231) are arbitrary the inequality £(t) >0 implies the positive semi-definiteness of the hermitian 
parts of the constitutive tensors in the Laplace domain: 


4> : t)'^ : 4>* > 0 

q p'^ q* >0 (25) 


where (p is an arbitrary complex-valued second order symmetric tensor, q is an arbitrary complex-valued vector, 
and the relation holds for all x and z (in the region of convergence). Similar results can be derived for the Fourier 
transform of the time derivatives of the constitutive tensors. Under the restriction that the support of a distribution 
f be bounded, its Fourier transform is given by (f (t), Now we let be equal to in Eq. (EiH) and 

follow the subsequent process to determine that the hermitian parts of the Fourier transforms of the time-derivative 
constitutive tensors must be positive semi-definite. However, the boundedness restrictions on the constitutive tensors 
need not be so severe for us to come to this conclusion. We merely assume that the constitutive tensors are distributions 
of slow growth to come to the same conclusion. To do so we consider the following for a test function 4){t) € S and a 
distribution f{t) G S : 

/ t poc -| pt pOO pOO 

ds(j)*{s) / dvf{v)(j){s — v) = — / ds(/)*(s) / da;/(a;) / (j){s — v)e~^‘^'^dv 

-oo J —oo J —oo J —oo J —oo 

^ J dujp{uj)f{uj)p*{uj) (26) 

where the distributional Fourier transform is used and the last step follows by choosing the test function 4>{t) such that 
it vanishes for t > t. It is clear from the above that under the much less restrictive conditions that the constitutive 
tensors be distributions of slow growth, Eq. o can be written, after some manipulations, in the following way: 


1 

27r 


/ OO pt 

d. 

-OO J —( 


ds(/)*(s)e / du(j)(u) 


Re 


du! 


r(w) : D(w) : cr*{uj) -|- u(a;) • p{uj) ■ u*(w) 


> 0 


(27) 


Since cr, u are arbitrary, the above will be satisfied only if the following is true about the hermitian parts of the Fourier 
transforms: 

0 ; D'* : 0* > 0 

q • • q* > 0 (28) 

Eqs. (I25I28I) constitute important conclusions in this paper and they are in line with earlier conclusions in circuit 
theory^. These relations can be used to place similar constraints on the transforms of constitutive tensors D, p. For 
this we need only consider the relations between the transforms of derivates as they apply to distributions in S . The 
following relations are noted: 


/ = -iw/; f = zf] f G s' (29) 

Since the elements of our constitutive tensors are assumed to be in S the above relations apply to them. Specifically 
we have, for instance for the Fourier transforms: 


D(a;) 


D(x,aj) 

iuj 


zD(x, Lu) 

UJ 


p(x,w) 


UJ 


(30) 
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Eq. (|30|) explicitly relates the Fourier transforms of the constitutive tensors to the Fourier transforms of their time 
derivatives. The Laplace transforms can be similarly related. To extend the constraints which passivity applies to 
the Fourier transforms of the constitutive tensors (analogous to Eq. 1281) we use the hermitian transpose operation, f. 
It is clear that and with similar relations holding for second order tensors. By expanding 

D, p, D, p into their hermitian and non-hermitial parts in Eq. (jSOp and by applying the hermitian transpose operator 
it becomes clear that the hermitian parts of the time derivative quantities are related to the non-hermitian parts of 
the original tensors in the following sense (w dependent implicit): 

D h ^-r^nh 'tji ^ ~nh \ 

= -D ; p = -p (31) 

i i 

Therefore, now Eq. (1^51) places the following constraints on the Fourier transforms of D, p: 

cj) : -D"''(x,a;) :cj)* >0 
i 

q-jP^’^M-q*>0 (32) 

It is clear that operations such as cj) : D"'^(x, uj) : cf)* and q-p”^(x, w) -q* result in purely imaginary numbers. However 
the factor i in the denominator ensures that the quantities in Eq. (1321) are purely real. Passivity and causality of the 
system demand that these numbers also be non-negative for positive values of uj (and non-positive for negative values 
of uj). We also note the corollary result that had we decided to represent Fourier transform through the exponential 
g-zcjt of 6*“^* we would have arrived at the complementary result where the non-hermitian quantities above 

would have been required to be negative semi-definite instead of positive semi-definite. In the following sections we 
will consider a specialization and a generalization of the above results. The specialization refers to cases where the 
constitutive tensors possess major symmetries and the generalization refers to the above results in the context of more 
general forms of constitutive relations such as the Willis kind of coupled relations. 


IV. WITH MAJOR SYMMETRIES 

Up to now we have assumed no special forms for the compliance and density tensors beyond the minor symmetries 
which ensure rotational stability of the system. We now consider the specialization of the above results to a case 
where the constitutive tensors possess major symmetries as well. For the density tensor we mean that its components 
satisfy = pji. Similarly we require that the fourth order compliance tensor iJyfei satisfy Dij^i = Since the 

components of the Fourier and Laplace transforms of the constitutive tensors are only related to the corresponding 
time domain components, it is clear that these major symmetries will extend to them as well. With these additional 

requirements Eqs. (IMl) imply that D^* = and p^* = p^ essentially meaning that and p^ are composed 
only of the real parts of D and p respectively. Similarly, and are composed only of the imaginary parts of 

D and p respectively. Furthermore, these characteristics should hold for all tensors of current interest (i.e. Fourier 
and Laplace transforms of D, p as well). Consideration of this specialization is of interest because for this case, the 
definiteness relations apply simply to the real and imaginary parts of the relevant tensors. Specifically, we have the 
following relations for this case: 

</> : Re D : </>* > 0; q • Re p • q* > 0 V</) = 0^, q 

cj) : —Im D :</>*> 0; q • —Im p • q* > 0 V</) = cp^, q (33) 

i i 

Restricting our attention for this section to the real parts (denoted by the subscript r) of the time derivative, fourier 

transformed tensors and to the imaginary parts (denoted by the subscript i) of the fourier transformed tensors, 

and using the shorthand > to imply positive semi-definiteness, the above relations are condensed to (uj dependence 
implied): 

Dr, ^ 0 V X, W 

Similar relations are derived by Milton and Willis^ in the context of minimum variational principles for time-harmonic 
waves in a dissipative medium. Since the imaginary parts of the Fourier transformed constitutive relations, for the 
simpler major-symmetric case as this one, corresponds to the dissipation in the system we note an interesting result as 
a corollary. For constitutive relations which do not necessarily possess major-symmetry, it is the non-hermitian parts 
of the Fourier transformed tensors which corresponds to dissipation in the system. In other words, a conservative 
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system can be expected to be hermitian in the Fourier transform of its constitutive relations. One system which 
immediately corresponds to the major-symmetric specialization being considered here is the case of one dimensional 
elastodynamics: 

€{x,t) = / ds D{x,t — s)a{x, s) p{x,t) = / ds p{x,t — s)u{x, s) 

J —oo —oo 

Causality and passivity results from the earlier sections can be immediately extended here to conclude that -D,/3 are 

analytic for Re z > 0 and that Dr, Pr, Dr, pr,ioDi, uipi > 0 where the symbol > actually means greater than or equal 
to in the present case and not just positive semi-definiteness. The static case, for which the constitutive tensors 
can be represented through the delta distribution, is seen to trivially satisfy these conditions since <5 = j = 1. This 
ensures that ujDi, ujpi = 0 etc. The symbol > can be understood to mean greater than equal to, and not just positive 
semi-definiteness, whenever a diagonal constitutive relation is being considered. In those cases passivity and causality 
would dictate that the > relations apply to the diagonal elements individually. Diagonal relations for the density 
tensor include those cases where ptj oc 6ij, and for the compliance tensor include those cases where Dijki oc SikSji. 

V. GENERALIZATION TO OTHER CONSTITUTIVE RELATIONSHIPS 

To derive the passivity relationships we required that energy could be expressed in a particular form (which it 
does automatically for the constitutive relations considered up to now). We will use this observation to generalize the 
results from the previous sections to more general constitutive relations such as the Willis relations. In the subsequent 
treatment we will understand the space dependence to be implicit in the sense of Eq. (HI. Let w(t),v(t) denote 
column vectors consisting of n time dependent tensors. Elements of w(t) are assumed to be in D to ensure that 
the energy expression exists (elements of v(t) are also distributions but they need not be so restricted). Let v(t) be 
derivable from w(x, t) through a linear, real, time invariant, and causal relationship v = L * w where L is a n x n 
matrix of real valued tensors: 


it) = f 

J —( 


ds 'Lit — s)w(s) = / ds L(s)w(t — s) 


(34) 


Each element of every tensor in L is assumed to be a distribution of slow growth. Let us also assume that the energy 
absorbed by the system up to a time t can be represented by: 

/ t pt poo 

ds w'l'(s)v(s) = Re / dsw'l'(s) / du L(u)u(s — u) (35) 

-OO J —OO J —oo 


then some conclusions apply to the Fourier and Laplace transforms of L and L. For instance, decomposing L into 
its hermitian and non-hermitian parts L = -|- L"^ and noting that (ylL^y)! = y’l'L^y and, therefore, real and 

(ytL^-'iyjt = —ytL"'*y and, therefore, imaginary. This means that £{t) emerges from (or L^). Passivity requires 
that absorbed energy must be non-negative at all times or that £{t) > 0. This implies the following: 

y^£v>0; y^LV>0; y^-L"'‘y>0 (36) 

I 

The Willis relations are a special case of the above. Specifically, they relate stress, strain, velocity, and momentum in 
a coupled constitutive relationship (x, t dependence implied): 



Power from Eq. (Ha can be written as: 

Pit) = Re [cTijeL -I-Pitt*] = Re [aljCij + u*Pi] = Re w’^'v 

with 



(37) 


(38) 


V = 


(39) 
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Since energy can now be written as 


£{t) = f ds P(s) = Re f ds w'l'(s)v(s) 
J —oo J —oo 


(40) 


we immediately have from the earlier results in this section that passivity implies that the following tensors will be 
positive semi-definite in the sense of Eq. (p6|) : 


L'* = 


D Si 

+ 

D Si 

t 

; = 

D Si 

+ 

1 

•>d-> 

, I/3-> 

t 

; L”'* = 

D Si 


D Si 

_^2 ~P_ 


>2 ~P_ 


k p. 


S2 p 


Sa p 


Sa p 


(41) 


It must be noted that the above results make no assumptions about the process by which Willis properties have been 
defined or derived or whether the relations are hermitian or not. They are merely the constraints which must be 
satisfied if an elastodynamic system, represented in the Willis form, is required to be causal. In fact, the question 
of whether the Willis constitutive relation displays self-adjointness (or are hermitian) or not has been addressed in 
literature. In Ref.'*^^ it has been shown that the property of self-adjointness is preserved at the level of the effective 
response. In other words, the effective relations are (not) self-adjoint if the constituting materials themselves are 
(not) self-adjoint. Ref4^ talks about the related question of hermiticity. The Willis relations exhibit several degrees 
of non-uniqueness^ and, to the author’s knowledge, it is not clear if some, many, or all sets of Willis properties (that 
can be assigned under any given case) satisfy the kinds of causality and passivity requirements discussed in this paper. 


VI. CONCLUSIONS 

In this paper we clarify the constraints that causality and passivity place on the elastodynamic constitutive tensors. 
Analogous questions have been addressed in other fields but the elastodynamic case is generally more complicated 
due to the higher order and non-diagonal nature of its constitutive relations. Here we deal with the problem in 
considerable generality wherein the elements of the constitutive tensors are assumed to be generalized functions in 
time. The treatment and conclusions presented here, therefore, apply to metamaterial applications which often involve 
singular and coupled constitutive forms and also to the static limit where the constitutive tensors are in the form of 
delta distributions. Specifically we show that the satisfaction of passivity (and causality) directly requires that the 
hermitian parts, as defined later, of the transforms (Fourier and Laplace) of the time derivatives of the elastodynamic 
constitutive tensors be positive semi-definite. Additionally, the conditions subsequently require that the non-hermitian 
parts of the Fourier transforms of the constitutive tensors be positive semi-definite for positive values of frequency 
and negative semi-definite for negative values of frequency. We show that when major symmetries are assumed these 
definiteness relations apply simply to the real and imaginary parts of the relevant tensors. For diagonal and one¬ 
dimensional problems, these positive semi-definiteness relationships reduce to simple inequality relations over the real 
and imaginary parts. Finally we extend the results to highly general forms of constitutive relations which include the 
Willis inhomogeneous relations as a special case. 
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